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r- H Abstract: In this paper, our focus is on certain applications for mobile robotic 

1 ''^ networks, where reconfiguration is driven by factors intrinsic to the network rather 

than changes in the external environment. In particular, we study a version of the 
coverage problem useful for surveillance applications, where the objective is to po- 

" sition the robots in order to minimize the average distance from a random point in 

^ a given environment to the closest robot. This problem has been well-studied for 

^\ omni-directional robots and it is shown that optimal configuration for the network is 

_> a centroidal Voronoi configuration and that the coverage cost belongs to ©(m,"^"), 

1"^ where m is the number of robots in the network. In this paper, we study this problem 

!/3 for more realistic models of robots, namely the double integrator (DI) model and 

I ^1 the differential drive (DD) model. We observe that the introduction of these motion 

constraints in the algorithm design problem gives rise to an interesting behavior. 

fSJ For a sparser network, the optimal algorithm for these models of robots mimics 

^ that for omni-directional robots. We propose novel algorithms whose performances 

OO are within a constant factor of the optimal asymptotically (i.e., as m -^ +00). In 

\l particular, we prove that the coverage cost for the DI and DD models of robots is 

'i-^ of order tti"^'^. Additionally, we show that, as the network grows, these novel algo- 

^ ; rithms outperform the conventional algorithm; hence necessitating a reconfiguration 

f in the network in order to maintain optimal quality of service. 
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1 Introduction 

The advent of large scale sensor and robotic networks has led to a surge of 
interest in reconfigurable networks. These systems are usually designed to 
reconfigure in a reactive way, i.e., as a response to changes in external condi- 
tions. Due to their importance in sensor network applications, reconfiguration 
algorithms have attracted a lot of attention, e.g., see [9]. However, there are 
very few instances in engineering systems, if any, that demonstrate an inter- 
nal reconfiguration in order to maintain a certain level of performance when 
certain intrinsic properties of the system are changed. However, examples 
of endogenous reconfiguration or phase transitions are abound in nature, e.g.. 
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desert locusts [4^ who switch between gregarious and social behavior abruptly, 
etc. An understanding of the phase transitions can not only provide insight 
into the reasons for transitions in naturally occurring systems but also iden- 
tify some design principles involving phase transition to maintain efficiency in 
engineered systems. 

In this paper, we observe such a phenomenon under a well-studied setting 
that is relevant for various surveillance applications. We consider a version 
of the so-called Dynamic Traveling Repairperson Problem, first proposed by 
[T^ and later developed in [5]. In this problem, service requests are generated 
dynamically. In order to fulfill a request, one of the vehicles needs to travel 
to its location. The objective is to design strategies for task assignment and 
motion planning of the robots that minimizes the average waiting time of a 
service request. In this paper, we consider a special case of this problem when 
service requests are generated sparingly. This problem, also known as coverage 
problem, has been well-studied in the robotics and operations research com- 
munity. However, we consider the problem in the context of realistic models 
of robots: double integrator models and differential drive robots. Some pre- 
liminary work on coverage for curvature-constrained vehicles was reported in 
our earlier work [7]. In this paper, we observe that when one takes into con- 
sideration the motion constraints of the robots, the optimal solution exhibits 
a phase transition that depends on the size of the network. 

The contributions of this paper are threefold. First, we identify an interest- 
ing characteristic of the solution to the coverage problem for double integrator 
and differential drive robots, where reconfiguration is necessitated intrinsically 
by the growth of the network, in order to maintain optimality. Second, we pro- 
pose novel approximation algorithms for double integrator robots as well as 
differential drive robots and prove that they are within a constant factor of 
the optimal in the asymptotic (m -^ -l-oo) case. Moreover, we prove that, 
asymptotically, these novel algorithms will outperform the conventional al- 
gorithms for omni-directional robots. Lastly, we show that the coverage for 
both, double integrator as well as differential drive robots scales as 1/m^^^ 
asymptotically. 



2 Problem Formulation and Preliminary Concepts 

In this section, we formulate the problem and present preliminary concepts. 

Problem Formulation 

The problem that we consider falls into the general category of the so called 
Dynamic Traveling Repairperson Problem, originally proposed in [121. Let Q C 
K^ be a convex, compact domain on the plane, with non-empty interior; we will 
refer to Q as the environment. For simplicity in presentation, we assume that 
Q is a square, although all the analysis presented in this paper carries through 
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for any convex and compact Q with non-empty interior in M^. Let A be the 
area of Q. A spatio-temporal Poisson process generates service requests with 
finite time intensity A > and uniform spatial density inside the environment. 
In this paper, we focus our attention on the case when A ^ 0+, i.e., when the 
service requests are generated very rarely. These service requests are identical 
and are to be fulfilled by a team of m robots. A service request is fulfilled 
when one of m robots moves to the target point associated with it. 

We will consider two robot models: the double integrator model and the 
differential drive model. The double integrator (DI) model describes the dy- 
namics of a robot with significant inertia. The configuration of the robot is 
g = {x,y,Vx,Vy) C M'' where {x,y) is the position of the robot in Cartesian 
coordinates, and (ua;,Wj,) is its velocity. The dynamics of the DI robot are 
given by 

x{t) = Vx{t), 

Vyit) == Uyit), U,x{tf + Uy{tf < U^^^^ Vt , 

where Wmax and Umax are the bounds on the speed and the acceleration of the 
robots. 

The differential drive model describes the kinematics of a robot with two 
independently actuated wheels, each a distance p from the center of the robot. 
The configuration of the robot is a directed point in the plane, g — {x, y, 9) C 
SE(2) where {x, y) is the position of the robot in Cartesian coordinates, and 
9 is the heading angle with respect to the x axis. The dynamics of the DD 
robot are given by 

X{t) = -{Wi{t)+ Writ)) cos e{t), 

y{t)^]^{wi{t) + Wr{t))sm9{t), 

9{t) = --{Wr{t) - Wi{t)), \wi{t)\ < Wmnxyt, \Wr{t)\ < WmaxVt, 

2p 

where the inputs wi and Wr are the angular velocities of the left and the right 
wheels, which we assume to be bounded by Wmax- Here, we have also assumed 
that the robot wheels have unit radius. 

The robots are assumed to be identical. The strategies of the robots in 
the presence and absence of service requests are governed by their motion 
coordination algorithm. A motion coordination algorithm is a function that 
determines the actions of each robot over time. For the time being, we will 
denote these functions as tt = (tti , 7r2, . . . , TTm), but do not explicitly state their 
domain; the output of these functions is a steering command for each vehicle. 
The objective is the design of motion coordination algorithms that allow the 
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robots to fulfill service requests efficiently. To formalize the notion of efficiency, 
let Tj be the time elapsed between the generation of the j-th service request, 
and the time it is fulfilled and let T^ := limj^+oo lim;s^^o+ E[Tj] be defined 
as the system time under policy tt, i.e., the expected time a service request 
must wait before being fulfilled, given that the robots follow the algorithm 
defined by tt. We shall also refer to the average system time as the coverage 
cost. Note that the system time T^ can be thought of as a measure of the 
quality of service collectively provided by the robots. 

In this paper, we wish to devise motion coordination algorithms that yield 
a quality of service (i.e., system time) achieving, or approximating, the the- 
oretical optimal performance given by Topt — inf ^r T^^ ■ Since finding the opti- 
mal algorithm maybe computationally intractable, we are also interested in 
designing computationally efficient algorithms that are within a constant fac- 
tor of the optimal, i.e., policies tt such that T,^ < nTopt for some constant 
K. Moreover, we are interested in studying the scaling of the performance of 
the algorithms with m, i.e., size of the network, other parameters remaining 
constant. Since, we keep A fixed, this is also equivalent to study the scahng 
of the performance with respect to the density m/ A of the network. 

We now describe how a solution to this problem gives rise to an endogenous 
reconfiguration in the robotic network. 

Endogenous Reconfiguration 

The focus of this paper is on endogenous reconfiguration, that is a reconfig- 
uration necessitated by the growth of the network (as the term 'endogenous' 
implies in biology), as opposed to any external stimulus. We formally describe 
its meaning in the context of this paper. In the course of the paper, we shall 
propose and analyze various algorithms for the coverage problem. In particu- 
lar, for each model of the robot, we will propose two algorithms, tti and 712- 
The policy tti closely resembles the omni-directional based policy, whereas the 
novel 7r2 policy optimizes the performance when the motion constraints of the 
robots start playing a significant role. We shall show that 

lim — ^ ~ 1, lim — ^ = 0, limsup — ^ < c for some constant c > 1. 

This shows that, for sparse networks, the omni-directional model based al- 
gorithm TTi is indeed a reasonable algorithm. However, as the network size 
increases, there is a phase transition, during which the motion constraints 
start becoming important and the 1:2 algorithm starts outperforming the tti 
algorithm. Moreover, the 1:2 algorithm performs with a constant factor of the 
optimal in the asymptotic (to — > -l-oo) case. Hence, in order to maintain ef- 
ficiency, one needs to switch away from the tti policy as the network grows. 
It is in this sense that we shall use the term endogenous reconfiguration to 
denote a switch in the optimal policy with the growth of the network. 
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3 Lower Bounds 

In this section, we derive lower bounds on the coverage cost for the robotic 
network that are independent of any motion coordination algorithm adopted 
by the robots. Our first lower bound is obtained by modeling the robots as 
equivalent omni-directional robots. Before stating the lower bounds formally, 
we need to briefly review a related problem from computational geometry 
which has direct consequences for the case omni-directional robots. 

The Continuous m-median Problem 

Given a convex, compact set Q C M^ and a set of points p — {pi,P2, ■ ■ ■ ,Pm} G 
Q™, the expected distance between a random point q, sampled from a uniform 
distribution over Q, and the closest point in p is given by 

Hm(p, Q):= / — min ||p, - g|| dq = V / —\\p^-q\\dq, 
jQAie{i,...,m} T^i-JVdp)'^ 

where V{p) — (Vi(p), V2(p), • ■ • , Vm(p)) is the Voronoi (Dirichlet) parti- 
tion [TU] of Q generated by the points in p, i.e., 

V»(p) = {qe Q: \\q~P^\\ < ||g-Pj||,Vje {!,..., m}}, ie {!,..., m}. 

The problem of choosing p to minimize Ti^ is known in geometric optimiza- 
tion [T] and facility location [6 literature as the (continuous) 77i-mcdian prob- 
lem. The TO-median of the set Q is the global minimizcr 

P*miQ) ^ argmiupgg^ Hmiv, Q)- 

We let Ti^{Q) — Hm(pJ^(Q), Q) be the global minimum of Tim- The solution 
of the continuous m-median problem is hard in the general case because the 
function p i—^ Timip, Q) is not convex for m > 1. However, gradient algorithms 
for the continuous multi-median problem can be designed [5]. We would not 
go further into the details of computing these 77i-median locations and assume 
that these locations are given or that a computationally efficient algorithm 
for obtaining them is available. 

This particular problem formulation, with demand generated indepen- 
dently and uniformly from a continuous set, is studied thoroughly in |llj for 
square regions and J13j for more general compact regions. It is shown in |13| 

that, in the asymptotic {m — > -|-cx)) case, T-L^{Q) = Chcx\/ almost surely, 
where Chcx ~ 0.377 is the flrst moment of a hexagon of unit area about its 
center. This optimal asymptotic value is achieved by placing the m points on 
a regular hexagonal network within Q (the honeycomb heuristic). Working 
towards the above result, it is also shown in [13 that for any to G N: 

•^ <Kn(Q)<c(Q)./?, (1) 

TTTO, V TO, 
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where c(Q) is a constant depending on the shape of Q. In particular, for a 
square Q, c(Q) « 0.38. 

We use two different superscripts on Topt for the two models of the robots, 
i.e., we use r5,j for the DI robot and T^^ for the DD robot. Finally, we state 
a lower bound on these quantities as follows. 

Lemma 1. The coverage cost satisfies the following lower bound. 

°P' - V ' °P* - w 

Proof. The proof follows trivially by relaxing the constraints on the robots 
and allowing them to move like omni-directional robots with speeds Wmax for 
DI robots and Wmax for DD robots. One can then adopt the lower bound on 
coverage cost for omni-directional robot, e.g., [2] to arrive at the result. 

Remark 1. Since from Equation (fTl, Ti.^{Q) G i7(l/\/m). Lemma fl] implies 
that T^i and T°d also belong to i7(l/^M)- 

This lower bound will be particularly useful for proving the optimality of 
algorithms for sparse networks. We now proceed towards deriving a tighter 
lower bound which will be relevant for dense networks. The reachable sets of 
the two models of robots will play a crucial role in deriving the new lower 
bound. We study them next. 



Reachable Sets for the Robots 

In this subsection, we state important properties of the reachable sets of the 
double integrator and differential drive robots that are useful in obtaining 
tighter lower bound. 

Let r : G X M^ ^ M+ be the minimum time required to steer a robot 
from initial configuration g in G to a point q in the plane. For the DI robot, 
G = K'* and for DD robots G — SE(2). With a slight abuse of terminology, 
we define the reachable set of a robot, TZt{g), to be the set of all points q in 
Q that are reachable in time t > starting at configuration g. Note that, in 
this definition, we do not put any other constraint (e.g., heading angle, etc.) 
on the terminal point q. Formally, the reachable set is defined as 

ntig) = {qeR''\Tig,q)<t}. 

We now state a series of useful properties of the reachable sets. 

Lemma 2 (Upper bound on the small-time reachable set area). 

The area of the reachable set for a DI robot starting at a configuration 
g — {x, y, i, y) G M.*, with i^ + y^ — Vq, satisfies the following upper bound. 

Area(7^,(5)) < ( ^^o^m^.i^ + o{t'), as t ^ 0+ for v, > 0, 
\<axt for Vo = 0. 
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The area of the reachable set for a DD robot starting at any configuration 
g € SE(2) satisfies the following upper bound. 

Area(7^t(g)) < ^"'L.^' + o(t'), as t ^ 0+. 

Proof. The result for the differential drive robot has been derived in [5]. We 
derive the result for the double integrator robot here. Assume, without any 
loss of generality, that the robot is initially placed at the origin with ve- 
locity aligned with the x-axis. The maximum of the absolute value of the 
x-coordinate and y-coordinate of all the points reachable in time t is less than 
or equal to v^t + ^Umaxi'^ and ^Umaxi^, respectively. Therefore, the area of 
the reachable set is trivially upper bounded by 4{vot + ^Wmax^^JCjWmaxi'^) = 

2woMmaxt^ + "niax^"'- 

Lemma 3 (Lower bound on the reachable set area). The area of the 
reachable set of a DI robot starting at a configuration g — {x,y,x,y) G M^, 
with x"^ + y'^ = vq, satisfies the following lower bound. 

Area(7^,(g)) > ^^t^ Vt < f ^. 

The area of the reachable set of a DD robot starting at any configuration 
g £ SE(2) satisfies the following lower bound. 

Area(7^t(5)) > ^^i^.t^. 

Lemma 4. The travel time to a point in the reachable set for a DI robot 
starting at a configuration g = (x, j/,i,y) G M*, with i^ + y^ = Vq, satisfies 
the following property. 



.^g^q)dq>'^^^^II^t^ Vt 



< 



TT V, 



TZtig) "Umax 

The travel time to a point in the reachable set for a DD robot starting at any 
configuration g € SE(2) satisfies the following property. 

[ r{g,q)dq>'^t\ 
Proof. For both the robots, /.j, , s t((?, q)dq — L Area,{TZs{g))ds . Using Lemma 3 

2 — 

for a double integrator robot, for all t < 5^^^ we have that, 

f / \J \ WQUmax /"' 3j t^0Wmax,4 

/ T{g,q)dq> / g-'ds^ ^——t . 

JTZtig) •^ "'0 ^^ 

The proof for the differential drive robot follows along similar lines. 
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We are now ready to state a new lower bound on the coverage cost. 

Theorem 1 (Asymptotic lower bound on the coverage cost). The cov- 
erage cost for a network of DI or DD robots satisfies the following asymptotic 
lower bound. 

liminf TZm'/' > ^(.^ ) '^\ and 



6pA\^/ 



m^+oo ^ OWrn.,^ V 5 / 



3 



Proof We state the proof for the double integrator robot. The proof for the 
differential drive robot follows along similar lines. 

In the following, we use the notation A; = Area(I?Vi(5)), where 'DVi{g) := 
{f? e Q I T{g,,q) < T{gj,q) Vj 7^ i}. We begin with 



-DI _ • r Y^ ' 1 



> inf / — min T(qi,q)dq. (2) 

Let RAiigi) be the reachable set starting at configuration g and whose area 
is Ai. Using the fact that, given an area A^, the region with the minimum 
integral of the travel time to the points in it is the reachable set of area A^, 
one can write Equation (pi) as 



m ^ 

- ^S! 1^ / ^^(ft,?) dq. (3) 



rpDI 

-'opt 



Let ti be defined such that Area(7?,t . (gi)) = A^. Lets assume that as m — > +cx), 
Ai -^ 0+ (this point will be justified later on). In that case, we know from 

n / A \^/^ • 

Lemma 2 that, ti can be lower bounded as ti > i ^a^u — ) ' where Vi is the 

speed associated with the state gi. Therefore, from Equation pi and Lemmal4J 
one can write that 

m 

^o°pt> M E/ -A9^,q)dq 

> inf yl^'^ll'll^tf. (4) 



Using the above mentioned lower bound on ti, Equation Bl) can be written as 
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mil A 



4/3 



11 " 

~24v/2yt,,^/3 ,//3 {Ai.A2,....A„}eR'"^^ * 

subject to 2_. Ai = -4 and Aj > Vi e {1, . . . , m}. 

Note that the function f{x) = x*''^ is continuous, strictly increasing and 
convex. Thus by using the Karush-Kuhn- Tucker conditions [3 , one can show 
that the quantity X^T-'^i ^^ minimized with an equitable partition, i.e., 
Ai = A/m, Vi. This also justifies the assumption that for m -^ +oo, A^ -^ 0+. 

Remark 2. Theorem fl] shows that both T^^^ and T^^ belong to n{l/m^/^). 

4 Algorithms and their Analyses 

In this section, we propose novel algorithms, analyze their performance and 
explain how the size of the network plays a role in selecting the right algorithm. 
We start with an algorithm which closely resembles the one for omni- 
directional robots. Before that, we first make a remark relevant for the analysis 
of all the algorithms to follow. 

Remark 3. Note that if rip is the number of outstanding service requests at 
initial time, then the time required to service all of them is finite (Q being 
bounded). During this time period, the probability of appearance of a new 
service requests appearing is zero (since we are dealing with the case when 
the rate of generation of targets A is arbitrarily close to zero). Hence, after an 
initial transient, with probability one, all the robots will be in their stationary 
locations at the appearance of the new target. Moreover, the probability of 
number of outstanding targets being more than one is also zero. Hence, in the 
analysis of the algorithms, without any loss of generality, we shall implicitly 
assume that there no outstanding service requests initially. 

The Median Stationing (MS) Algorithm 

Place the m robots at rest at the m-median locations of Q. In case of the 
DD robot, its heading is chosen arbitrarily. These ?n-median locations will be 
referred to as the stationary locations of the robots. When a service request 
appears, it is assigned to the robot whose stationary location is closest to 
the location of the service request. In order to travel to the service location, 
the robots use the fastest path with no terminal constraints at the service 
location. In absence of outstanding service tasks, the robot returns to its 
stationary location. The stationary configurations are depicted in Figure [l] 
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Fig. 1. Depiction of typical stationary configurations for the Median Stationing 
Algorithm. On the left: Dl robots at rest at their stationary locations. On the right: 
DD robots with arbitrary headings at their stationary locations. In both the figures, 
the shaded cell represents a typical region of responsibility for a robot. 



Let r^g and T^ be the coverage cost as given by the above policy for 
the DI and DD robots, respectively. 

Theorem 2 (Analysis of the MS algorithm). The coverage cost for a 
network of DI and DD robots with the MS algorithm satisfies the following 
hounds. 






I2V2A 



2u,, 



nUQ) ^ rr^BD ^ K.(S) 



^ rriUU ^ 



piT 



2Wr. 



Proof. The lower bounds on T^ and T^g follow trivially from Lemma n] 

For a double integrator robot, the minimum travel time from rest at loca- 
tion p to a point q is given by 



r((p,0),g)< 



nv-q\\ 



\\p-q\\ 



for ||p-g|| < ^ 
otherwise. 



Therefore, for any g, r((p, 0), q) can be upper bounded as 



^((P,0),<7)< 



\\P-<1\\ 



inp-4 



< 



\\p-q\\ 



I2V2A 
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The upper bound on T^g is then obtained by taking the expected value over 
all g e Q while taking into consideration the assignment policies for the 
services to robots. For a differential drive robot, the travel time from any 
initial configuration {p, 9) to a point q can be upper bounded by ^~ + ^w'' ■ 
The result follows by taking expected value of the travel time over all points 
inQ. 

Remark 4- Theorem^ and Lemma IT] along with Equation ([T| imply that, 

lini ±Mf = l and lim ^ = 1. 

This implies that the MS algorithm is indeed a reasonable algorithm for sparse 
networks, where the travel time for a robot to reach a service location is almost 
the same as that for an omni-directional vehicle. 

However, as the density of robots increases, the assigned service locations 
to the robots start getting relatively closer. In that case, the motion con- 
straints start having a significant effect on the travel time of the robots and 
it is not obvious in that case that the MS algorithm is indeed the best one. 

In fact, for dense networks, one can get a tighter lower bound on the 
performance of the MS algorithm: 

Theorem 3. The coverage cost of a dense network of DI and DD robots with 
the MS algorithm satisfies the following bounds: 



limmiT^lm^/'^ > 0.321 ' 

m — *oo 



liminfrDD>^^ 



1/4 



Proof. Let us consider the DI case first. The minimum travel time for the i-th 
robot to reach a point q inside its region of responsibility Vi, starting from 
rest at the median location pi is bounded by 

T{{pi,0),q)>maxi—===, \, (5) 



where 



dv; = max \\pi — q\\ 



For large m, it is known that the honeycomb heuristic is optimal |13| , yielding 



/ 2j^ 
lim dv,m^/2 = J^ w0.62\/^, Vi G {1, . . . , m}. (6) 



rn—> oo 



Clearly, in these conditions the first term in (pi is the dominant one. 
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Another consequence of the optimaHty of the honeycomb heuristic is that 
hm lril,{QW'^ ^ ci,,^VA. (7) 

m — >oo 

Integrating ([5| over Q, multiplying both sides by m}'\ and taking the limit 
as 77T, — > oo, we get 



liminf TDim'/" > Hminf . ^"'^^^"^ 



1/2 



•-MS 



"--"o V2M,„axrfv.ml/2' 



Finally, using ^ and ([t]), we get the desired result. 

We will only sketch the proof for the DD case. The minimum travel time 
for a DD robot can be decomposed into the sum of the cost of turning towards 
the target point, plus the Euclidean distance between the robot and the target 
point. The Euclidean term vanishes as m increases. The turning cost on the 
other hand remains bounded away from zero. Since the robot's initial heading 
is chosen randomly, the expected turning angle is 7r/4, which combined with 
the maximum turning rate ui^ax/p yields the stated result. 

In other words, for DI robots, the MS algorithm requires them to stay sta- 
tionary in absence of any outstanding service requests. Once a service request 
is assigned to a robot, the amount of time spend in attaining the maximum 
speed iimax becomes significant as the location of assigned service requests 
start getting closer. Similar arguments hold for DD robots. 

An alternate approach, as proposed in the next algorithm, is to keep the 
robots moving rather than waiting in absence of outstanding service requests. 
The algorithm assigns dynamic regions of responsibility to the robots. 

The Strip Loitering (SL) Algorithm 

This algorithm is an adaptation of a similar algorithm proposed in [7J for 
Dubins vehicle, i.e., vehicles constrained to move forward with constant speeds 
along paths of bounded curvature. 

Let the robots move with constant speed v* = min{w,„ax, 3 22""" i ^^'^ 
follow a loitering path which is constructed as follows. Divide Q into strips 

of width w where w = min ( f :^ A+w.ZBp'■^/A \ 2p*\, where p* := ^^ . 

Orient the strips along a side of Q. Construct a closed path which can be 
traversed by a double integrator robot while always moving with constant 
speed V* . This closed path runs along the longitudinal bisector of each strip, 
visiting all strips from top-to-bottom, making U-turns between strips at the 
edges of Q, and finally returning to the initial configuration. The m robots 
loiter on this path, equally spaced, in terms of path length. A depiction of the 
Strip Loitering algorithm can be viewed in Figure [2J Moreover, in Figure [3] 
we define two distances that are important in the analysis of this algorithm. 
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Variable d2 is the length of the shortest path departing from the loitering path 
and ending at the target (a circular arc of radius p*). The robot responsible for 
visiting the target is the one closest in terms of loitering path length (variable 
di) to the point of departure, at the time of target-arrival. Note that there 
may be robots closer to the target in terms of the actual distance. However, 
we find that the assignment strategy described above lends itself to tractable 
analysis. 
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Fig. 2. Depiction of the loitering path for the double integrator robots. The segment 
providing closure of the loitering path (returning the robots from the end of the last 
strip to the beginning of the first strip) is not shown here for clarity of the drawing. 



point. 

or 

departure 



Fig. 3. Close-up of the loitering path with construction of the point of departure 
and the distances J, di, and d-2. for a given target, at the instant of appearance. 



After a robot has serviced a target, it must return to its place in the 
loitering path. We now describe a method to accomplish this task through 
the example shown in Fig. [3] After making a left turn of length ^2 to service 
the target, the robot makes a right turn of length 2d2 followed by another left 
turn of length di, returning it to the loitering path. However, the robot has 
fallen behind in the loitering pattern. To rectify this, as it nears the end of 
the current strip, it takes its U-turn early. 
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Let Tg^ be the coverage cost as given by the above algorithm for DI robots. 
We now state an upper bound on T^j^. 

Theorem 4 (Analysis of the SL algorithm). The coverage cost for a team 
of DI robots implementing the SL algorithm satisfies the following asymptotic 
upper bound. 

lim sup TgD^mi/^ < ^{p*A+l0.38p*'VAy^\ 

m — >+oo V 

Proof Since a similar algorithm for Dubins vehicle was analyzed in [7], we 
only outline the proof here and refer to [7] for more details. Denote the length 
of the closed path as Li. Due to equal spacing of the robots along the loitering 
path, 

EMi]-|^. (8) 

We now calculate an upper bound on Li. To that effect, let A'strips be the 
number of strips, Lstrip be the length travelled along a single strip, Lu-turn 
be the length of a u-turn and Lciosure be the length of the closure path. With 
these notations, ii can be bounded as 

^1 — -^ *strips-^strip ' \-^^stiips -^J-^u— turn ~r -^closure- wj 

One can compute bounds for various terms on the right side of Equation d9| . 
Substituting these bounds into Equation (l9| and taking into account Equa- 
tion (Is]) we get that 

^^^^^^A+10Mp*VA^2VA + 6.19p*^ (10) 

2mw m 

To calculate Fi[d2] we define S as the smallest distance from the target to any 

point on the loitering path (see Fig. 

s < p* and 6 is uniformly distributed between and w/2 



31. Since ^2(5) = 2p*sm-\J^) for 



EM2] = i^ r 'sin-i (,r^) ds < ^y^. (11) 



w Jo VV 2p* J ~ 4 

Therefore, the coverage cost is given by 

Ts. < ^M±m\, (12) 

V* 

Therefore, from Equations (12 1, (10 1 and (ITT]) we get that for w < 2p, 

^SL < ^ i + i + y^Vp*w. (13) 

In order to get the least upper bound, we minimize the right hand side of 
Equation (13 1 with respect to w subject to the constraint that w < 2p*. 
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Remark 5. Theorem 4 and Theorem llimply that TS belongs to 0{l/m^^^). 
Moreover, Theorem 5 and TheoremTS] together with Equation IT] imply that 
^sl/^ms — > 0+ as to ^- +00. Hence, asymptotically, the SL algorithm out- 
performs the MS algorithm and is within a constant factor of the optimal. 



We now present the second algorithm for DD robots. 



The Median Clustering (MC) Algorithm 



Form as many teams of robots with k :- 



4.09 



(A) 



2/3 



,1/3 



robots in each 



team. If there are additional robots, group them into one of these teams. Let 
n :— [^J denote the total number of teams formed. Position these n teams at 
the n-median locations of Q, i.e., all the robots in a team are co-located at the 
median location of its team. Within each team j, j € {1, . . . , n}, the headings 
of the robots belonging to that team are selected as follows. Let ij > fc be the 
number of robots in team j. Pick a direction randomly. The heading of one 
robot is aligned with this direction. The heading of the second robot is selected 
to be along a line making an angle j, in the counter-clockwise direction, with 
the first robot. The headings of the remaining robots are selected similarly 
along directions making ^-angle with the previous one (see Figure l4| . These 
headings will be called the median headings of the robots. Each robot in 
a team is assigned a dominance region which is the region formed by the 
intersection of double cone making half angle of ^ with its median heading 
and the Voronoi cell belonging to the team (see Figure B . When a service 
request appears, it is assigned to the robot whose dominance region contains 
its location. The assigned robot travels to the service location in the fastest 
possible way and, upon the completion of the service, returns to the median 
location of its team and aligns itself with its original median heading. 






Fig. 4. Depiction of the Median Clustering Algorithm with teams of 4 robots each. 
The shaded region represents a typical region of responsibility for a robot. 
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Let T^c ^^ ^^^ coverage cost as given by the above policy. We now state 
an upper bound on T^^ in the following theorem. 

Theorem 5 (Analysis of the MC algorithm). The coverage cost for a 
team of DD robots while implementing the MC algorithm satisfies the following 
asymptotic upper bound. 

limsupT°gmi/3 < ^^{pA^/^ 

?Tl— > + CXD Wyjiax 

Proof. The travel time for any robot from its median location p to the location 
g of a service request is upper bounded by ^' + 2w'^ k ■ Taking the expected 
value of this quantity while taking into consideration the assignment policy 
of the service requests gives us that 

r-<M2) + ^^. (14) 

From Equation (fTl), 7i* (Q) < 0.38w ^. Moreover, for large m, n w m/k. This 



combined with Equation (14 1, one can write that, for large m, 



,DD ^ 0.38 I A , TTp 



^^ Wmax V m/k 2Wmaxfc ' 



,2/3 



The right hand side of Equation ( 15 1 is minimized when k = 4.09 1 -^ J 
Substituting this into Equation (151, one arrives at the result. 

Remark 6. Theorem 5 and Theorem 1 imply that T^^ belongs to 0{\/m}'^). 
Moreover, Theorem S and Theoreni|3] together with Equation [T] imply that 
-^McZ-^MS* ^ ^^ ^^ "^ ^ +00. Hence, asymptotically, the MC algorithm 
outperforms the MS algorithm and is within a constant factor of the optimal. 

5 Conclusion 

In this paper, we considered a coverage problem for a mobile robotic network 
modeled as double integrators and differential drives. We observe that the op- 
timal algorithm for omni-directional robots is a reasonable solution for sparse 
networks of double integrator or differential drive robots. However, these al- 
gorithms do not perform well for large networks because they don't take into 
consideration the effect of motion constraints. We propose novel algorithms 
that are within a constant factor of the optimal for the DI as well as DD 
robots and prove that the coverage cost for both of these robots is of the 
order 1/m}/^. 

In future, we would like to obtain sharper bounds on the coverage cost 
so that we can make meaningful predictions of the onset of reconfiguration 
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in terms of system parameters. It would be interesting to study the prob- 
lem for non-uniform distribution of targets and for higher intensity of arrival. 
Also, this research opens up possibilities of reconfiguration due to other con- 
straints, like sensors (isotropic versus anisotropic), type of service requests 
(distributable versus in-distributable), etc. Lastly, we plan to apply this re- 
search in understanding phase transition in naturally occurring systems, e.g., 
desert locusts [4]. 
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